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A significant observational effort has been directed to unveil the nature of the so-called dark 
energy. However, given the large number of theoretical possibilities, it is possible that such a task 
cannot be performed on the basis only of the observational data. In this article we discuss some 
thermodynamic properties of this energy component assuming a general time-dependent equation- 
of-state parameter oj(z) = wo + u a f{z), where uiq and ui a are constants and /(z) may assume 
different forms. We show that very restrictive bounds can be placed on the wq - w a space when 
current observational data are combined with the thermodynamic constraints derived. 

PACS numbers: 98.80.-k, 95.36.+X, 95.30.Tg 



I. INTRODUCTION 

The observational evidence for an acceleration of the 
expansion of the Universe is now overwhelming, although 
the precise cause of this phenomenon is still unknown 
(see, e.g., [l[ for recent reviews). In this concern, and be- 
sides the need for more accurate estimates of cosmological 
parameters, the current state of affairs also brings to light 
some other important aspects regarding the physics of the 
mechanism behind cosmic acceleration. Certainly, one of 
these aspects concerns the thermodynamical behavior of 
a dark energy-dominated universe, and questions such as 
"what is the thermodynamic behavior of the dark energy 
in an expanding universe?" or, more precisely, "what is 
its temperature evolution law?" must be answered in the 
context of this new conceptual set up. Another interest- 
ing aspect in this discussion is whether thermodynamics 
in the accelerating universe can place constraints on the 
time evolution of the dark energy and can also reveal 
some physical properties of this energy component. 

The aim of this paper is twofold. First, to derive phys- 
ical constraints on the dark energy from the second law 
of thermodynamics and to deduce the temperature evolu- 
tion law for a dark component with a general equation-of- 
state (EoS) parameter w(a); Second, to perform a joint 
statistical analysis involving current observational data 
together with the thermodynamic bounds on w(a). To 
do that, we assume the following generalized formula for 
the time evolution of w(a) [2| 



w(a) 



Px 
Px 



wq + w a 



Wq + W a 



I -a' 3 



P 

!-(! + *)' 



(1) 



which recovers some well-known EoS paramcterizations 
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in the following limits: 
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where p x and p x stand for the dark energy pressure and 
energy density, respectively (see also [6| for other EoS pa- 
rameterizations) . The analyses are performed using one 
the most recent type la supernovae (SNe la) observa- 
tions, the nearby + SDSS + ESSENCE + SNLS + Hub- 
ble Space Telescope (HST) set of 288 SNe la discussed 
m Ref. (which we refer to as SDSS compilation). We 
consider two sub-samples of this latter compilation that 
use SALT2 Q and MLCS2k2 SN la light-curve fit- 
ting method. Along with the SNe la data, and to help 
break the degeneracy between the dark energy parame- 
ters we also use the baryonic acoustic oscillation (BAO) 
peak at zbao = 0.35 [10 ] and the current estimate of the 
CMB shift parameter K = 1.71 ± 0.019 [lj. We work in 
units where c = 1. Throughout this paper a subscript 
stands for present-day quantities and a dot denotes time 
derivative. 



II. THERMODYNAMIC ANALYSIS 

Let us first consider a homogeneous, isotropic, spatially 
flat cosmologies described by the Friedmann-Robertson- 
Walker (FRW) flat line element, ds 2 = dt 2 - a 2 (t)(dx 2 + 
dy 2 + dz 2 ). where a(t) = 1/(1 + z) is the cosmological 
scalar factor. The matter content is assumed to be com- 
posed of baryons, cold dark matter and a dark energy 
component. 

In such a background, the thermodynamic states of a 
relativistic fluid are characterized by an energy momen- 
tum tensor (perfect-type fluid) 



a particle current, 
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and an entropy current, 

S a = nau° 



(2c) 



where h a @ := g a/3 — u a u^ is the usual projector onto 
the local rest space of u a and n and a are the particle 
number density and the specific entropy (per particle), 
respectively jT^j. The conservation laws for energy and 
particle number densities read 



'iff 



Px + (px +Px)9 = 



N a ; Q = ri + n6 = 



(3a) 
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where semi-colons mean covariant derivative, 9 = 3a/ a is 
the scalar of expansion and the quantities p x , p x , n and 
a are related to the temperature T trough the Gibbs law: 
nTda — dp x — Px + Pj: dn. From the energy conservation 
equation above, it follows that the energy density for a 
general w(a) component can be written as 



p x oc exp 
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da 



(4) 



Following standard lines (see, e.g., [T3, EI), ^ i s P os ~ 
sible to show that the temperature evolution law is given 
by 
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where we have split the dark energy pressure as 
Px = Pa + II = w p x + w a f(a)p x , 



(5) 



(6) 



where /(a) = (1 — ar)//3. By combining the above equa- 
tions, we also find that 



T oc exp 



w(a) 



da 



(7) 



where we have used that n oc a , as given by the con- 
servation of particle number density [Eq. ([3b| ]. For pa- 
rameterization P^, shown in Eq. (TTJ), Eq. ([7]) can be 
rewritten as 



T x = T£a 



-3(i»o+ 



exp 



3w a 
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(8) 



which reduces to the generalized Stefan-Boltzmann law 
for w a = [H (see also pjj). Fr om the above expres- 
sions, we confirm the results of Ref. [15| (for a constant 
EoS parameter) and find that dark energy becomes hot- 
ter in the course of the cosmological expansion since the 
its EoS parameter must be a negative quantity. A pos- 
sible physical explanation for this behavior is that ther- 
modynamic work is being done on the system (see, e.g., 
Fig. 1 of Ref. 15]). In particular, for the vacuum state 
(w = —1) we obtain T oc a 3 . To illustrate this behavior, 
Fig. 1 shows the dark energy temperature as a function 
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FIG. 1: Temperatute evolution law for parameterizations P2 
(fi = 0) and P3 (fi = 1) assuming some arbitrary values of 
Wo, and T§ = 10- 2 T CMB , where T CMB = 2.73 K. For 
comparison, we also show the CMB temperature curve (solid 
line). 



of the scale factor for f3 = (P2) and /3 = 1 (P3) 1 by as- 
suming arbitrary values of w , w a and T X = 1§~ 2 Tq MB , 
where Tq b — 2.73 K. From this analysis, it is clear 
that an important point for the thermodynamic fate of 
the universe is to know how long the dark energy temper- 
ature will take to become the dominant temperature of 
the universe. A basic difficulty in estimating such a time 
interval, however, is that the present-day dark energy 
temperature has not been measured, being completely 
unknown. 

By considering that the chemical potential for this 
u>(a)-fluid is null (as occurs for w = 1/3), the Euler's 
relation defines its specific entropy, i.e., 



N 



P., 



nT 



(9) 



Now, by combining the above equations, it is straight 
forward to show that the product p x o? /T = const 
that 



so 



S x oc (1 + w) 



(10) 



For a constant EoS parameter, the above expression re- 
covers some of the results of Ref. 15]. Note also that the 
vacuum entropy is zero (w — —1) whereas for phantom 
dark energy (w < 1), which violates all the energy con- 
ditions [17[, the entropy assumes negative values being, 



1 As is well known, PI has the drawback of becoming rather un- 
physical for too > 0, with the dark energy density p x blowing up 
as e 3w ° z at high-z. We, therefore, do not consider this parame- 
terization in our analyses. 
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FIG. 2: Contours of A\ 2 = 6.17 in the parametric plane wo - 
w a for P2 (/? = 0) . The light gray area represents the thermo- 
dynamic constraint of Eq. f) whereas the small hachured 
area is the resulting parametric space when the constraint 
from the second law of thermodynamics [Eq. (|12[) ] is added 
to the analysis. 



therefore, meaningless (For a discussion on the behav- 
ior of a phanton fluid with nonzero chemical potential, 
see [l8[ . See also [lj| for an alternative explanation in 
which the temperature of the phantom component takes 
negative values and 20] for other thermodynamic analy- 
ses of dark energy). 

Two cases of interest arise directly from Eq. (|10l) . The 
case in which S x — const, implies necessarily that wp — 
for all the above parameterizations 2 . The second case is 
more interesting, with the w(a)-fluid mimicking a fluid 
with bulk viscosity where the viscosity term is identified 
with the varying part of the dark energy pressure II. For 
this latter case, we note that the positiveness of S x im- 
plies that 



w a > 



1 + w 
/(<*) 



(11) 



which clearly is not defined at a — 1, where w = wq. 
Finally, by combining Eqs. © and (IT0|) with the conser- 
vation of the particle number density shown earlier, we 
obtain from the second law of thermodynamics that 



S a ; a <xw > 
or, equivalently, w a < 0. 



(12) 



2 Note that time-dependent EoS parameterizations without the 
constant term uiq are imcompatiblc with the case S x = constant. 



FIG. 3: The same as in Figure 2 for parameterization P3 
08 = 1). 



III. CONSTRAINTS ON THE w - w B PLANE 

In this section, we combine the above physical con- 
straints (fTT|) and (|12[) with current observational data in 
order to impose bounds on the dark energy parameters. 
We use one of the most recent SNe la data sets avail- 
able, namely, the SDSS compilation discussed in Ref. 0] ■ 
This compilation comprises 288 SNe la and uses both 
SALT2 |f and MLCS2k2 [<| light-curve fitters (see also 
21] for a discussion on these light-curve fitters) and is 
distributed in redshift interval 0.02 < z < 1.55. Along 
with the SNe la data, and to help break the degeneracy 
between the dark energy parameters wq and wp we use 
the BAO [13] and shift parameters [Til ] 



A = D\ 



H 2 



ZBAO 



= 0.469 ±0.017 , 



K = 0^f( ZC M B ) - 1.71 ± 0.019 , 



(13a) 



(13b) 



where Dy — [^(zbacO-zbaoZ-HX-zbao)] 1 ^ 3 is the so-called 
dilation scale, defined in terms of the dimensionless co- 
moving distance r, zbao = 0.35 and zcmb = 1089. In our 
analyses, we minimize the function \ 2 = Xsno + Xbao + 
XcmB' which takes into account all the data sets men- 
tioned above and marginalize over the present values of 
the matter density f2 m and Hubble parameters Hq. 

Figures 2 and 3 show the main results of our joint anal- 
yses. We plot contours of A% 2 = 6.17 in the parametric 
space wo - w a for P2 and P3, respectively. The light 
gray region displayed in the plots stands for the physi- 
cal constraint (|11[) . Since this inequality is a function of 
time, the region is plotted by assuring its validity from 
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a = 10~ 4 up to today at a = 1. The resulting para- 
metric space, when all the observational data discussed 
above are combined with the constraints ([TTj) and (fl~2"|) . 
corresponds to the small hachured area right below the 
w a = line. These results clearly illustrate the effect 
that the thermodynamic bounds discussed in the previ- 
ous section may have on the determination of the dark 
energy EoS parameters. In particular, we note that the 
resulting allowed regions are even tighter for the loga- 
rithmic parameterization P2 than for P3 (CPL). Since 
the SALT2 compilation allows for more negative values 
of wq, the joint constraints involving this SNe la sub- 
sample are also more restrictive (Figs. 2b and 3b). For 
completeness, we display in Table I the changes in the 
2a estimates of wq and w a due to the thermodynamic 
bounds HI]) and (fT2jl . 



IV. CONCLUSIONS 



In spite of its fundamental importance for an actual un- 
derstanding of the evolution of the universe, the relevant 
physical properties of the dominant dark energy compo- 
nent remain completely unknown. In this paper we have 
investigated some thermodynamic aspects of this energy 
component assuming that its constituents are massless 
quanta with a general time-dependent EoS parameter 
w(a). We have discussed its temperature evolution law 
and derived constraints from the second law of thermo- 
dynamics on the values of wo and w a for a family of w(a) 
parameterizations given by Eq. ([TJ. When combined 
with current data from SNe la, BAO and CMB observa- 
tions, we have shown that such constraints provide very 
restrictive limits on the parametric space wq - w a (see 
Figs. 2 and 3). 

Finally, it is also worth mentioning that in the present 
analysis we have assumed that the chemical potential fi 
for the w(a)-fluid representing the dark energy is null. A 
more general analysis relaxing this condition (p =/= 0) is 
currently under preparation and will appear in a forth- 
coming communication. 
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